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Abstract 

In this paper we discuss the evolution of domain walls generated in the 
early Universe considering, unlike the previous studies, an interaction 
between the walls and a major gaseous component of the dark matter. 
The walls are supposed able to reflect the particles elastically and with 
a reflection coefficient of unity. We discuss a toy Lagrangian that could 
give rise to such a phenomenon. In the simple model studied we obtain 
highly non-relativistic and slowly varying speeds for the domain walls 
( ^ 10~ 2 (1 + z) ~ 1 ) and negligible distortions of the microwave background. 
In addition, these topological defects may provide a mechanism of form- 
ing the large scale structure of the Universe, by creating fluctuations in 
the dark matter Spj p ~ 0(1) on a scale comparable with the distance 
the walls move from the formation ( in our model d < 20 h 1 Mpc). The 
characteristic scale of the wall separation can be easily chosen to be of 
the order of 100 Mpc instead of being restricted to the horizon scale, as 
usually obtained. 
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1. Introduction 




The cosmological consequences of primordial phase transitions associated with 
scalar fields have been the subject of many studies in recent years. The topologi- 
cal defects created in the transitions, such as domain walls, strings and monopoles 
are potentially of great interest for Cosmology, since they could supply seeds for the 
formation of the large scale structure of the Universe. Specifically, domain walls are 
sheet-like regions of false vacuum in-between domains having different and discon- 
nected vacuum ground states of the scalar field. The simplest and the most studied 
model involves a real scalar field with a quartic potential and a negative sign for 
the mass term; after the phase transition the field rolls down to one of the two zero 
temperature minima for the potential; this leaves a domain structure on scales bigger 
than the correlation length of the field, resembling closely what happens in an Ising 
model 1 . When originally introduced, the phase transitions considered were on the 
GUT scale 2 . The trouble is that domain walls on the GUT scale rapidly become the 
dominant form of matter in the Universe and produce much too big distortions in the 
present microwave background. 

Recently, interest in domain walls has been raised again considering late phase 
transitions (at z ~ 100) that would give rise to so called ’’soft’’ domain walls 3 . These 
walls may never be massive enough as to distort the microwave background but may 
a priori be a dominant gravitational component of the present Universe, triggering 
the formation of galaxies and changing the expansion rate. These possibilities have 
been excluded by a numerical study 4 of the evolution of the field itself through the 
phase transition and after, as the walls appear and evolve by their surface tension . 
The domain walls soon reach relativistic speeds and the average scale of the system 
becomes comparable to the horizon scale, making these walls unusable for the forma- 
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tion of the large scale structures we "see 5 ' 6 . Very similar results have been obtained 7 
by considering directly the evolution of the walls after the phase transition . In that 
calculation the approximation taken is that the wall thickness is much smaller than 
the radius of curvature of the wall surface. 

The problems mentioned arise due to the lack of energy dissipation in the models 
considered; the mass-energy stored in the walls gets efficiently converted into their 
kinetic energy, rapidly raising them to relativistic speeds. We therefore consider the 
effect of introducing in the equation of motion of the walls a friction term that is a 
function of the wall speed relative to the background matter and to its density. The 
idea of studying the consequences of friction on domain walls can be traced back to 
refs. 1,8 ,9, but it was never fully developed because it was introduced in the context of 
GUT scale phase transitions, in which case inc lu ding friction would even worsen the 
problems pointed out previously. In this paper we will consider much lower energy 
scales, of the same order of those obtained in ref.3. It will be shown that indeed 
there exists an interesting range of the wall energy density for which the average 
” inter- wall ’’distance is of the order of 100 Mpc today, and that these domain walls 
are compatible with the limits on the anisotropy of the microwave background. 

The paper is organized in the following way: in section 2 we derive the equation 
of motion of an element of domain wall without any friction term other than the 
usual due to the universal expansion; in section 3 we concentrate our attention on 
the friction pressure arising when walls move through a homogeneous gas reflecting 
all incident particles elastically; in section 4 we introduce the results of section 3 into 
the equation of motion previously calculated; in section 5 we discuss what kind of 
particle Lagrangian may lead to the premises of this paper and the consequences of 
our model on the microwave background 
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2. The equation of motion in the absence of friction 


In order to approach the problem we will assume, first of all, that we are deal- 
ing with domain walls late enough after the phase transition so that the thin wall 
approximation can be considered roughly valid 7 . We are therefore interested in the 
motion of sharp interfaces under their surface tension. The shape of the network, 
which is related to the details of the model chosen, will turn out to be unimportant 
in the discussion. The important assumption made here is that we are dealing with 
two or more degenerate values for the vacuum ground state, so that the driving force 
for the motion of the walls is only due to their surface tension. 

Many approaches may be considered to get the local equation of motion of the 
walls; the most direct of these is just to start with the well known equation of motion 
for real scalar fields 4 ’ 10 , 


a- 1 -72,- dV 


$ + 3-$ = - 




(1) 


where a is the scale factor, given by a = f 2/3 if fi = 1 ( we express fin 2tf 0 J /3 units, 
H 0 being the present value of Hubble constant; a = 1 today). Eq. (1) is expressed 
in comoving coordinates and universal time. After the phase transition there are 
regions of different vacua separated by kinks (which are classical solutions of eq.(l)). 
Throughout the following calculations we will assume that these kinks are moving 
non-relativistically; this will turn out to be a sensible choice ( see section 4 ). 

In general we can define as the surface of a kink the 2-d space on which dV/d$ = 0; 
at each point of the surface we can label as x the axis normal to it. If x 0 (t) is the 
intercept of the surface with the x axis and the principal radii of curvature at the 
point are much bigger than the wall thickness A we can represent the kink by a 
function $(x - x 0 (f)). The calculation is easily performed when we recall that 
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where x is the unit vector perpendicular to the kink surface. The divergence V -x 
is, in 3-d, the sum of the curvatures along the principal axes of the surface at the 
point considered n . In this way from eq. (1) we get 
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Evaluating (2) at x = x 0 we see that d$jdx becomes very big when A — ► 0 
{d$/dx ~ A^ 1 ), while the d 2 $jdx 2 term is very small ( it would be exactly 
d 2 $ /dx 2 = 0 if the wall were straight ). In the thin wall approximation we therefore 
get to the final expression 

1 '±+i 

‘.**1 r 2 . 

where rj and r 2 are the principal comoving radii of curvature at the chosen point of 
the wall network. 

If there were no universal expansion (set a = 1 constant ) then eq.(3) would 
look like x 0 = — (1/rj + l/r 2 ) ; if a is the mass-energy density of the walls, then 
P T = a(\/ri -j- l/r 2 ) would simply be the pressure due to the surface tension, exactly 
the same form that one obtains in condensed matter. This also reminds us that eq.(3) 
is just Newton’s second law divided by a. 


3. The friction term 


We now derive the pressure exerted on the walls moving with speed 1 through 
some homogeneous medium interacting with it. We are going to study only the case 
in which the medium remains homogeneous throughout all the period of evolution 
considered. This can be considered valid, for example, if perfectly reflecting walls 


move so little that they are not able to reshuffle the bulk of the matter, i.e. if they 
move of a small fraction of the distance between each other, so that no particle 
interacts with two different walls in a cosmological time (see also the discussion at 
page 7). In all the following we are restricting ourself to this simple case. 

We begin by writing down the general expression for the friction force acting on 
the domain walls as they move non-relativistically through a homogeneous gas. The 
particles coupled with the walls will be taken to be weakly interacting (WIMP’s). We 
will also suppose that the walls are able to reflect elastically all the incident particles, 
regardless of their energy at the impact. This condition could be relaxed, as we will 
discuss at the end of this section. 

For a non-relativistic gas we can write that the pressure exerted by the gas on the 
wall is given by (see Appendix): 

roc 

Pj - - 2 mn I B~ 2a {y-y x ) f(\y x \) dy x + 

2 mn T B- 2a {y-y x j* f{\y x \)dy x (4) 

J — oo 

where B = m/7 1 , y = B a v and y x = B a v x . 

Let’s consider the limits in which y <C 1 and y 1 . In the former case the 
thermal speed of the particles is much greater than the speed of the domain wall, 
since the average thermal momentum of the particles is p ~ T\ in the latter case the 
wall is moving through particles effectively at rest and the volume spanned remains 
depleted of the gas. The case y » 1 will turn out to be the most interesting in our 
discussion. 

For y <C 1, changing the variables inside the integrals ( y\ = y — yx] and expanding 
f(\y x \) in power series around j/ 1? we get 

Pj = —4mnB~ 2a y y\f\y\)dyi = 4 mn{T/m) Q vF (5) 
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where F = — / 0 °° y\ f'{y\) dy x is a constant of order unity. For y » 1 instead we 
obtain 


Pj = 2 mn v 2 (6) 

since f{y x ) ~ 0 for y x ;» 1 and therefore we can substitute (y — y x ) 2 by y 2 in the 
integrals. Such a result is not surprising if we recall that in this case A p ~ 2mv and 
that the number of particles hitting the unit area per unit time is nv. 

The case in which the gas is relativistic is even easier, since the number of particles 
hitting the wall per unit time is simply given by n ( c = 1 ) on both sides of the surface. 
Taking v <C 1 we get 


Pf = n 


J r oo f oo 

l Px(l +v) f{\p x \)dp x - / Px(l - v)f(\p x \)dp x 
0 JO 


= 2nvp x 


(?) 


so that Pj ^ vT 4 which is the limit discussed in the review 1 . 

W'e can also give an evaluation of the average thermal speed of the particles in- 
teracting with the walls e.g. for light neutrinos and for gas having the Boltzmann 
distribution, supposing that the particles decouple at a certain z d (note that through- 
out the paper 2 + 1 = a"* 1 ). The momentum of the particles shifts with the expansion 
of the Universe so that p d hmv(z) = (z d + 1 )/(z + 1), where 7 ” 3 = y/l - F 2 ; if at 
decoupling T d > m then p d ~ T d \ if T d <£ m then mv 2 d ^ T d , so that p d ~ (mT) 1/2 , 
Assuming the particles to be non-relativistic today we get 
f T d \'t 2 -z + 1 


for Td <C m and 


Zd + 1 


(8) 


u(?) 


i) 


z + 1 
Zd + 1 


(9) 


for Td » m. If we assume neutrinos of mass m — 10 eV and Td ~ 1 Mev we get 
F(ir) ~ 10 -5 at z = 0. This result will turn out to be useful in the following discussion. 
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In closing this section we go back briefly to our initial assumptions. Although 
we are interested here in studying the consequences of a reflection coefficient close to 
unity, there could be cases in which one has to deal with an energy dependent partial 
transmission of the incident particles through the walls. This would lead to a class of 
solutions in which the walls may decouple from the matter after a certain stage, when 
their speed with respect to the matter becomes bigger than a critical value. These 
possibilities are at present under investigation and go beyond the goals of this paper. 


4. Domain walls and friction: a simple case 


In this section we introduce the previously calculated friction term into the equa- 
tion of motion of the domain walls. 

Define R = ax and v = ax Q . We can rewrite eq.(3) in terms of physical coordinates 
and peculiar speed v and add up the friction term we have been talking about in the 
previous section. As we are going to show soon the most interesting case to study is 
the y > 1, when the domain walls move fast with respect to the thermal motion of 
the gas. If the matter interacting is a major component of the dark matter then it 
follows that the gas is highly non relativistic and y > 1 during most of the evolution 
of the network, as shown in a check of self-consistency at the end of this section. 

We define p m = mn so that p m = p mo /a 3 , where p mo is the mass density of the 
matter interacting with the walls today. Inserting eq.(6) into eq.(3) we get 12 


„ a ki 2 
v -I- 2-v + — v 2 = 


( 10 ) 


where k\ = 2p mo /<7. The constant a is the energy density of the walls. Let’s now 
suppose that v < k^v 2 /a 3 and 2 a/a < hv/a 3 . This also will turn out to be self- 
consistent. We therefore remain with the important equation: 
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w— K {-k + &) 


( 11 ) 


where we define K = k l i . 

Eq.(ll) can be finally written in the more useful comoving coordinates as 



( 12 ) 


We are interested in how much the walls move in average from the original con- 
figuration, and therefore we average the curvature over the surface S of the network 
contained in an arbitrarily large volume V r 3 . The quantity f is the mean interwall 
distance in comoving coordinates, defined as a point by point average of the distance 
to the next neighbouring wall integrated 13 over S. We therefore obtain 




,ri r 2 . 

In general r = 0r (with 0 ~ 0(1) ), so that we finally get to 

(5J) 1/2 = 


(13) 


K_ 

Qf 


1/2 


(14) 


where K — a/2p mo . 

Since we are studying the case in which f changes little from the phase transition 
to the present time, we can consider eq.(14) as an estimate of the average comoving 
speed r, which remains roughly constant. Our goal would be to determine a given r, 
but before doing so we should slightly modify eq.(14) taking into account the following 
correction. The friction term we utilize in eq.(10) is based on the assumption that 
only one reflection occurs to each particle. In comoving coordinates the speed of a 
free particle goes like r p ~ a ~ 2 . In the case y >> 1 after being reflected particles 
have a speed double than that of the wall, but this decreases due to the universal 
expansion and soon they get scattered again. A priori this fact could modify the form 
of the law of motion of the walls, but this is not the case, as we will see in a moment. 
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Let’s give a numerical estimate. Consider a wall moving at constant comoving speed 
r (we will check the validity of the assumption at the end of the calculation). At 
a certain time A a particle at rest is reflected and its speed is henceforth given by 
r p (t) = 2 r(L/<) 4/3 (since a = t 2/3 if fi = 1). The maximum comoving distance x max 
from the wall is reached at the time f mai , when the wall and the particle have equal 
speed. This yields r p (t ma x) = r —> t max = 1.7 1{ — * x max — 0.2rt max = 0.2ra ^ , which 
is roughly as far as any scattered particle can get from the kink. 

All the matter in the volume swept by the wall from its formation is contained 
within a distance x max in front of the kink, while the total distance traveled by the 
wall is Ar ~ rt max = ra 3/2 . Since A r/x max ~ const., we take the density of the 
matter in front of the wall to be roughly constant; in this way we get pj ron t ~ 6 p m at 
any time. This means that the initial assumption of constant comoving speed Eq.(14) 


is, at least approximately, self-consistent, substituting A by I\' = A/6, 

/ „ V ' 2 

(ii)' n = 


(15) 


A more accurate calculation would require a numerical simulation that takes into 
account, also gravitational effects. 

To proceed w r e need now an estimate of A''. If we associate the walls with the peaks 


of the distribution of galaxies observed in the survey 6 , which suggests the domain 
walls may be related to the clustering process, then the scale of our network will be 

, — n 

f = R 0 = 120 h~ l Mpc today. In this way we obtain K' = o/\2p m0 = 6/3 ■ 10 -2 Ar ; 
this can also be written as 

~ 12 /3~Kr 2 (16) 

since Q wo /^mo ~ cr/poRo by geometry. If we want the walls to produce a density 
fluctuation on a scale of the order, for example, of 20 A -1 Mpc in a major component 
of the dark matter (assuming Q mo = 1), setting 'Sr = 10 -2 = 20 h~ l Mpc yields 
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n wo ~ 1.2/9 10 “ 3 (equivalently, a ~ 1.2 0MeV 3 ); this constitutes an upper bound on 
n wo in our model and shows once again that the domain walls never get to dominate 
the energy density of the Universe. 

We can now easily see that the self-consistency conditions on eq.(12) are ensured. 
In fact, we know that 


(v 2 ) 1 / 2 ~ alO 


-2 


(H) 


which says that our initial assumption v <C 1 is satisfied by a big margin. It also 
means that at z > 30 the speed of the walls is smaller than the average thermal speed 
of particles like neutrinos with mass m > 10 eV (see estimate given by eq.(9) and 
compare it with eq.(17)). In general, when this occurs the friction term we used eq.(6) 
has to be replaced by eq.(5). At early times one should write (/J mo /cro 3 ) ( T 0 /am ) v ~ 
l jar — > v ~ 20 a 3 = 20 (2 + l) -3 (where T 0 = T to j ay ) but, since the evolution of the 
configuration takes place for the most part at 2 < 5, such a change wouldn’t affect 
our previous conclusions. 

Let’s continue our self-consistency check, going back to our original eq.(10). We 
know that in a Universe with critical density Q mo = 1 and a — t 2/3 , so that vjv = 2/3 1 
at each point of the network; we therefore get 


10 4 


,2 


u=--<12(^^ N )i = 

3 1 V <7 RJ 2 6/9 • 10 -2 1 2 


v > 4 • 10" 6 ^ 


and 


h 1.5 -10V 

2 ~ v < ^77 * v > 8 • 10 ~ 6 (3t 

a pt* 


As anticipated at the beginning of this section, at any time considered we meet 
the conditions for friction dominated motion (analogous considerations apply also at 
z > 30 in the m = 10 eV neutrino case). 
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5. Discussion 


We have shown that, if domain walls coming from some late phase transition are 
able to perfectly reflect the particles of gas of a component of the dark matter, then 
the domain wall network is bound to expand with the scale factor, provided that 
Q, wo /Q, mo < 1.2/310 -3 . The coupling between the scalar field $ and the particles in 
question (call the associated field $) may assume the very simple form of a mass 
term dependent on the spatial coordinates. For the sake of discussion take 'P to be 
fermions. A toy Lagrangian for 'P could be written as = 'P ^'P + (m + /($))'P'P 

(a Lagrangian of this form is obtained e.g. in ref.3 where /($) is a real function of the 
field 4> that gives the domain structure and it is assumed to get higher values within 
the kinks than outside. Due to the presence of the kinks, the mass of the particles 
changes when they get close to the solilon. In the non-relativistic case this situation 
is equivalent to obtaining a Schroedinger equation for free particles with a potential 
V = /(<P(r)). We can say that V(f) is a perfect barrier if the reflection coefficient on 
both sides of the kink is unity. 

Now we turn to consider a possible estimate of the wall thickness, that has been, 
up to this point, a free parameter. If the domain walls maintain their position from the 
formation (in the comoving coordinate system), we are actually bound to consider 
second order phase transitions not earlier than zj = R 0 /Rh{ z j ) = Ro/3if — 2 
10 ~ 2 z 3 / 2 -+ Zj < 2500 ( R h is the horizon scale at z } ) due to simple causality- 
considerations. As a consequence there is a lower bound to the thickness A of the 
domain walls; since the interwall distance is Rj ~ A ~ R 0 zJ 1 at formation, A > 
3.5 • 10 -5 in our units, which is A > 7 • 10 -2 /i -1 Mpc. Such a distance is far greater 
than the wavelength usually associated to any dark matter particle candidate ( e.g. 
for neutrinos A thermal < 10 4 d /_1 at any z). We infer that F(r) can be considered a 
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classical barrier of height E max such that for E < E max the reflection coefficient is 
unity and for E > E max it is zero. In this paper we only consider the case E max — ► oc. 

A couple of issues still remain to be solved. The walls carry a gravitational field 
that shifts the frequency of the microwave background radiation a slight amount w T hen 
this passes through the potential. Such a problem has been treated in refs. 3,14 . 

The infinitesimal shift of the average photon energy T while the photon is moving 
for a dt time through the gravitational influence of a wall is given by dT 6T + TSV , 

where V is the gravitational potential of the wall. V is roughly given by V ~ GoR 
at a distance R from the kink surface, within a cut off value ~ R/2; R is the average 
interwall distance in physical coordinates at the time considered 1S . The value of V 
varies in time due to the evolution of the network, so that 8V = ( dV/dt)6t + VV • 

^ -f 

<5/?(where |<57?| = 6t). We want to calculate the total shift in the temperature of the 
photons as they pass through the gravitational potential of a single wall, i.e. within 
the cut off distance of V . If we take roughly dV/dt ~ GcR in a region of order R in 
size(this is clearly an over-estimate), when we integrate the above expression for dT 
to find the total shift of the temperature we get a term 8T/T = ctGoR 2 /t in addition 
to the usual term due to the expansion; q is a fudge factor of order unity and t is the 
age of the Universe at the epoch considered. The biggest distortion can be reached 
at the present epoch: 8T/T ~ 10 -8 . 

Another effect may be considered. The fluctuation in the matter density due to 
the sweeping action of the wall gives rise to a gravitational influence limited to the 
region of thickness d ~ 20h~ 1 Mpc in which 6p m /p m f 0. The minimum value of the 
gravitational potential just due to this distribution of matter is V' m ~ Gp m (P . Using 
the same arguments as above we can calculate the distortion due to the matter in 
6T/T\ m ~ fdGpmod 3 jt {0 is fudge factor of order unity). Again the biggest 6T/T is 
reached today: 6TjT ~ 10 -7 [d/20/i -1 A/>c] 3 . 
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All other effects, including gravitational distortion at the last photon scattering 
surface (if zj > 1000 — ► 8T/T\iss ~ G<yR 0 ciLSSi with o.lss = 10 ) and effects 

originated at the phase transition, are comparatively much smaller. 

The values obtained refer to the distortion originated from a single wall. Even 
supposing that the phase transition takes place before the photon decoupling there 
are only N ~ Rh\R 0 ~ 3/6 • 10~ 2 = 50 walls between us and the last surface of 
scattering. An evaluation of the ST/T\ m due to the matter swept from the walls, 
which is the biggest distortion, can be obtained multiplying the single wall distortion 
by y/N and gives 6T/T\ m 10 6 . For the effects directly related to the domain 
walls our values of 6T/T are, for the same a, one order of magnitude lower than that 
calculated in the previous papers 8T/T ~ 10 aGoTPjt ~ 10 -7 ; this derives from an 
interwall separation an order of magnitude smaller. 

The gravitational interaction of the domain walls with matter is secondary with 
respect to the sweeping action. In fact, taking for the sake of discussion the favorable 
case of straight infinite walls, the peculiar speed gained by the particles due to the 
gravitational influence would be, after a cosmological time 14 , of the order v m ~ 

2xG(7f ~ lO -45 . 

In concluding the discussion we point out that one can also consider late first 
order phase transitions in order to achieve our big values for the average interwall 
distance, even while starting with a much smaller comoving correlation length at the 
critical temperature. In this way one can remove the lower bound on A obtained in 
this section. Such an analysis is left for future investigation. 
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6. Conclusion 


This paper wants to offer a framework for future work. We have made the following 
assumptions: 

• A network of domain walls is established in the primordial Universe through a 
second order phase transition. 

• The walls interact with an important gaseous component of the present energy 
density of the Universe, reflecting elastically all incoming particles regardless of 
their kinetic energy. 

• The configuration is bound to expand with the background comoving coordi- 
nates (up to higher order corrections). 

We reach the following conclusions: 

• There is a wide range of values for a , surface density of the domain walls, such 
that R 0) the average inter- wall distance today is of the order of the large scale 
structure observed for the galaxies. 

• The mechanism that generates the fluctuations in the distribution of the dark 
matter could be related also to the particle Lagrangian, and not just gravita- 
tional. 

• This suggests that the large scale structure could indeed form in intimate con- 
nection with the presence of the domain walls, although studying the evolution 
of the fluct uati ons an d the long distance gravitational effects (see the discussion 
on the Great Attractor in ref.14 obtained goes beyond the present work. 

• Domain walls never come to dominate the energy density of the Universe. 
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• Walls with a of Mev order and such a small interwall separation ( R 0 ~ 100 Mpc) 
are not able to distort the microwave background. Also the effects related to 
the matter density fluctuations are small. 

Some of the assumptions made to obtain our results may be relaxed, giving rise to 
the different scenarios we earlier mentioned. Particularly intriguing is the possibility 
of the wall decoupling mentioned in section 3 16 : domain walls may give rise to a 
spectrum of density perturbations and at some point decouple and start growing in 
the way described in the previous work 4 . This paper represents just a first attempt 
to approach the late phase transition issue from an angle that could solve some of the 
problems other investigations have found. This paper is meant to stimulate interest 
in such non-standard scenarios. 
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Appendix 

Consider an infinite wall moving along the x axis of a chosen coordinate system, 
with speed v > 0. At one side of the wall particles with speed i? x get reflected gaining 
momentum A p = 2m (v — v x ) (particles with v x > v will not interact with the wall). 
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The momentum distribution of the particles is defined on a 3 + 3 dimensional phase 
space; neverthless we are interested in the statistical distribution of the momenta only 
in the x direction, and we therefore integrate out all other degrees of freedom. In 
this way we can write, in a very general way, a statistical distribution f((m/T) a \v x \) 
defined so that B° /( B° \v x |) dv x = 1 ( where B = m/T) u . The coefficient 

q depends on the actual original distribution we are considering. For a Boltzmann 
distribution a = 1/2 while for light neutrinos (m -C 1 Mev ) a = 1. 

There are dN = B a n{v — v x )f(B Q \v x \)dv x interactions per each second and per 
unit area with momentum exchange A p ( n is the number density of the particles). 
On the other side of the wall Ap has opposite sign, so that we can write that the 
pressure exerted by the gas on the wall is given by: 

too 

Pj = -2 mn B- 2a (y-y x ) 2 f(\y x })dy x + 

J y 

2 mn f B~ 2q (y - y x ) 7 f(\y r \)dy x 

J — CO 


where y = B Q v and y x = 5 a n x .The first integral refers to particles having speed 
v x > v > 0 and hitting the wall from the back, while the second refers to particles 
hitting the wall from the front. 

We now derive, as an example, the form that f{B Q v) assumes in the case of light 
neutrinos ( m <C IMeV). We start up with the statistical distribution of neutrinos 
in thermal equilibrium: 


" (2 1) 3 / 


exp 


d 3 p 

Wm 2 


+ 1 


at T > Tj. At T < Trf this becomes 

<p P 


"=( 2 WJ 


Tj + 1 


expy Jp- + 

where T = T d a/a d \ since m(T d < 1 we get 


10 


= 9 f —P 

(27r) 3 J exp(p/T ) + 1 
at all times. 

The probability of finding a particle in an interval p x ,p x + dp x of the x component 
of the momentum is then 

< \ ' 9 r 27r Px^i 

9 Px ~n{2ir)*Jo exp ^^K^ + 1 

where px is an Y component of the momentum perpendicular to p x . Changing 
variables we get 


f(V ) = _L_r (is) 

4 £( 3 ) Jo exp\Jy\ + j/x + 1 

w T hich is an implicit function of y x = mv x /T (j/x = P±/T). A similar calculation can 
be performed for a Boltzmann distribution. 
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